Quadrupole Collective Dynamics from Energy Density Functionals: Collective 
Hamiltonian and the Interacting Boson Model 



O 

(N 

C 
3 



K. Nomura, 1 ' 2 T. Niksic, 2 T. Otsuka, 1 ' 3 ' 4 N. Shimizu, 1 ' 3 and D. Vretenar 1 - 2 

1 Department of physics, University of Tokyo, Hongo, Bunkyo-ku, Tokyo, 113-0033, Japan 
2 Physics Department, Faculty of Science, University of Zagreb, 10000 Zagreb, Croatia 
3 Center for Nuclear Study, University of Tokyo, Hongo, Bunkyo-ku Tokyo, 113-0033, Japan 
4 National Superconducting Cyclotron Laboratory, Michigan State University, East Lansing, MI 

(Dated: January 20, 2013) 

Microscopic energy density functionals (EDF) have become a standard tool for nuclear structure 
calculations, providing an accurate global description of nuclear ground states and collective exci- 
tations. For spectroscopic applications this framework has to be extended to account for collective 
correlations related to restoration of symmetries broken by the static mean field, and for fluctua- 
tions of collective variables. In this work we compare two approaches to five-dimensional quadrupole 
dynamics: the collective Hamiltonian for quadrupole vibrations and rotations, and the Interacting 
Boson Model. The two models are compared in a study of the evolution of non-axial shapes in Pt 
isotopes. Starting from the binding energy surfaces of 192 > 194 ' 196 pt i calculated with a microscopic 
energy density functional, we analyze the resulting low-energy collective spectra obtained from the 
collective Hamiltonian, and the corresponding IBM-2 Hamiltonian. The calculated excitation spec- 
tra and transition probabilities for the ground-state bands and the 7- vibration bands are compared 
to the corresponding sequences of experimental states. 

PACS numbers: 21.10.Re,21.60.Ev,21.60.Fw,21.60. Jz 
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I. INTRODUCTION 

One of the major research topics in the theoretical nu- 
clear structure physics has been the study of quadrupole 
collective dynamics from a microscopic viewpoint [lHj]- 
Quadrupole collectivity results from multi-nucleon dy- 
namics of nuclear surface deformation. The equilibrium 
shape of a nucleus can change depending on the num- 
ber of valence nucleons: a spherical vibrator, a deformed 
rotor, or a soft shape in between. In most isotopic or iso- 
tonic sequences the transition between different shapes is 
gradual, but in a number of cases, with the addition or 
subtraction of only few nucleons, one finds signatures of 
abrupt changes in observables that characterize equilib- 
rium shapes. These structure phenomena have been in- 
vestigated using concepts of quantum shape/phase tran- 
sitions [H, and advanced self-consistent beyond- mean- 
field approaches d-Q. 

Microscopic studies based on energy density function- 
als (EDFs) have been quite successful in reproducing 
with remarkable accuracy various intrinsic (bulk) prop- 
erties of medium-mass and heavy nuclei such as binding 
energies, density distributions, charge radii, giant res- 
onances, etc 0, The current generation of EDFs 
includes non-relativistic Skyrme- [Tu [ill ] and Gogny- 
[l2l [l3j functionals, as well as relativistic density func- 
tionals [H], [H| • The framework of EDFs has also been 
extended beyond the mean-field level to describe excita- 
tion spectra and electromagnetic transition rates. Models 
have been developed that perform restoration of symme- 
tries broken by the static nuclear mean field, and take 
into account quadrupole fluctuations: GCM configura- 
tion mixing calculations 0, 1 EMI, 

and solutions of 

the collective Hamiltonian with quadrupole degrees of 



freedom 0,0, [HHI. 

Another successful approach to the low-lying structure 
of medium-mass and heavy nuclei is based on the inter- 
acting boson approximation |22j. The Interacting Bo- 
son Model (IBM), in particular, provides not only an 
algebraic but in some cases also a microscopic descrip- 
tion of nearly-spherical and 7- unstable shapes |22h27| . 
From a microscopic point of view, the collective J = + 
and 2 + pairs of valence protons (neutrons) are mapped 
onto the corresponding boson images with J — + and 
2+ denoted by (s„) and {d v ) bosons, respectively 
[23j . A number of studies have been carried out to derive 
the IBM Hamiltonian starting from nucleonic degrees of 
freedom in terms of the conventional mapping method 
that starts from the shell model [23[, the more recent 
approach based on binding energy maps calculated with 
microscopic EDFs [HI • etc. 

A static self-consistent mean-field solution in the in- 
trinsic frame, for instance a map of the energy surface 
as a function of quadrupole deformation, is characterized 
by symmetry breaking: translational, rotational, particle 
number, and can only provide an approximate descrip- 
tion of bulk ground-state properties. To calculate ex- 
citation spectra and electromagnetic transition rates in 
individual nuclei, it is necessary to include correlations 
that arise from symmetry restoration and fluctuations 
around the mean-field minimum. Both types of correla- 
tions can be included simultaneously by mixing angular- 
momentum projected states corresponding to different 
quadrupole moments. The most effective approach for 
configuration mixing calculations is the generator coor- 
dinate method (GCM), with multipole moments used as 
coordinates that generate the intrinsic wave functions. 
It must be noted that, while GCM configuration mix- 
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ing of axially symmetric states has been implemented by 
several groups and routinely used in nuclear structure 
studies [28l - l30l |. the application of this method to triax- 
ial shapes presents a much more involved and technically 
difficult problem [IH, HJ ■ In addition, the use of general 
EDFs, that is, with an arbitrary dependence on nucleon 
densities, in GCM type calculations, often leads to dis- 
continuities or even diver gen ces of the energy kernels as 
a function of deformation [32ll33j. Only for certain types 
of density dependence a regularization method can be 
implemented, which corrects energy kernels and removes 
the discontinuities and divergences (34T - [36l | . 

In an approximation to the full GCM approach to 
five-dimensional quadrupole dynamics that restores ro- 
tational symmetry and allows for fluctuations around 
the triaxial mean-field minima, a collective Hamilto- 
nian can be formulated, with deformation-dependent 
parameters determined by constrained microscopic self- 
consistent mean-field calculations. The dynamics of the 
five-dimensional Hamiltonian for quadrupole vibrational 
and rotational degrees of freedom is governed by the 
seven functions of the intrinsic quadrupole deformations: 
the collective potential, three vibrational mass parame- 
ters, and three moments of inertia for rotations around 
the principal axes [1, [t| [l9l - l2"lj . 

Another approximation consists in mapping the self- 
consistent mean-field solution to a boson (IBM) Hamil- 
tonian. In Refs. |26l.l27| the energy surface for quadrupole 
degrees of freedom, calculated from a microscopic EDF, 
was mapped onto the corresponding boson energy surface 
under certain approximations. The interaction strengths 
of the boson Hamiltonian are determined by the mapping 
procedure. One then proceeds to calculate the excita- 
tion spectra and wave functions in the laboratory frame 
[26l The validity of the method of Ref. (2(| was 

tested in various mass regions [23, [33, [38| ■ 

It would be, therefore, interesting to compare the two 
approximations starting from the same self-consistent 
mean-field solution based on a microscopic EDF. In this 
work we compare spectroscopic observables calculated 
with the IBM Hamiltonian to the solution of the col- 
lective quadrupole Hamiltonian, with both calculations 
based on relativistic Hartree-Bogoliubov (RHB) [1J] self- 
consistent binding energy surfaces. The framework of 
relativistic EDFs and the corresponding collective Hamil- 
tonian have successfully been employed in studies of the 
evolution of ground-state shapes and spectroscop ic prop - 
erties of medium- heavy and heavy nuclei fl. foL Ha. [2ll 139 ] . 
In the present analysis we consider the even-even isotopes 
192- 196 j n tfi e framework these 7-soft nuclei 

can be characterized by the 0(6) dynamical symmetry 

The paper is organized as follows. In Sec. [Til we briefly 
describe the theoretical procedures used to derive the col- 
lective Hamiltonian and the IBM Hamiltonian starting 
from a given EDF. The microscopic RHB energy sur- 
face and the mapped IBM energy surface are discussed 
in Sec. IIIIl Spectroscopic properties of 192 ~ 196 pt cal- 



culated with the two models are compared in Sec. IIVI 
Section [V] summarizes the results and presents a short 
outline of future work. 

II. THEORETICAL FRAMEWORK 

The map of the energy surface as a function of the 
quadrupole collective variables f3 and 7 [I] is obtained 
from self-consistent RHB calculations with additional 
constraints on the axial and triaxial mass quadrupole 
moments. The quadrupole moments can be related to 
the polar deformation parameters /3 and 7. The param- 
eter /3 is simply proportional to the intrinsic quadrupole 
moment, and the angular variable 7 specifies the type 
and orientation of the shape. The limit 7 = cor- 
responds to axial prolate shapes, whereas the shape is 
oblate for 7 = 7r/3. Triaxial shapes are associated with 
intermediate values 0<7<7r/3. In this work the con- 
strained RHB calculations have been performed using 
the relativistic functional DD-PCI [42| . Starting from 
microscopic nucleon self-energies in nuclear matter, and 
empirical global properties of the nuclear matter equa- 
tion of state, the coupling parameters of DD-PCI have 
been determined in a careful comparison of the calcu- 
lated binding energies with data, for a set of 64 axially 
deformed nuclei in the mass regions A sa 150 — 180 and 
A sa 230 - 250. DD-PCI has been further tested in a 
series of calculations of properties of spherical and de- 
formed medium- heavy and heavy nuclei, including bind- 
ing energies, charge radii, deformation parameters, neu- 
tron skin thickness, and excitation energies of giant mul- 
tipole resonances. For the examples presented here, pair- 
ing correlations have been taken into account by em- 
ploying a pairing force that is separable in momentum 
space, and is completely determined by two parameters 
adjusted to reproduce the empirical bell-shaped pairing 
gap in symmetric nuclear matter (43[. 

The entire dynamics of the collective Hamiltonian is 
governed by seven functions of the intrinsic deformations 
/? and 7: the collective potential, the three mass parame- 
ters: Bpp, -B/3 7 , i? 77 , and the three moments of inertia 2*. 
These functions are determined by the choice of a partic- 
ular microscopic nuclear energy density functional and a 
pairing functional. The quasiparticle wave functions and 
energies, that correspond to constrained self-consistent 
solutions of the RHB model, provide the microscopic in- 
put for the parameters of the collective Hamiltonian [2l[ : 

H co \\ = T vi b + T ro t + VcoU , (1) 



with the vibrational kinetic energy: 
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and rotational kinetic energy: 

fc=i K 

V co \\ is the collective potential. Jk denotes the com- 
ponents of the angular momentum in the body-fixed 
frame of a nucleus, and the mass parameters Bpp, Bp-,, 
B 77 , as well as the moments of inertia Ik, depend on 
the quadrupole deformation variables j3 and 7: Xfc = 
4_Bfc/3 2 sin 2 (7 — 2kn/3). Two additional quantities that 
appear in the expression for the vibrational energy: r — 
B 1 B 2 B^, and w = BppB 11 — B„ , determine the volume 
element in the collective space. The moments of iner- 
tia are computed using the Inglis-Belyaev (IB) formula 
pil l45j . and the mass parameters associated with the 
two quadrupole collective coordinates qo — (Q20) and 
12 = (Q22) are calculated in the cranking approxima- 
tion. The potential V co n in the collective Hamiltonian 
Eq. (fTJ) is obtained by subtracting the zero-point energy 
corrections from the total energy that corresponds to the 
solution of constrained RHB equations, at each point on 
the triaxial deformation plane. The Hamiltonian Eq. (TTJ) 
describes quadrupole vibrations, rotations, and the cou- 
pling of these collective modes. The corresponding eigen- 
value problem is solved using an expansion of eigenfunc- 
tions in terms of a complete set of basis functions that de- 
pend on the deformation variables /3 and 7, and the Euler 
angles (/>, 9 and i/> [HI]. The diagonalization of the Hamil- 
tonian yields the excitation energies and collective wave 
functions for each value of the total angular momentum 
and parity, that are used to calculate observables. An im- 
portant advantage of using the collective model based on 
self-consistent mean- field single- (quasi)particle solutions 
is the fact that physical observables, such as transition 
probabilities and spectroscopic quadrupole moments, are 
calculated in the full configuration space and there is no 
need for effective charges. Using the bare value of the 
proton charge in the electric quadrupole operator, the 
transition probabilities between eigenvectors of the col- 
lective Hamiltonian can be directly compared with data. 

In an equivalent approach the RHB binding energy sur- 
face can be mapped onto the IBM Hamiltonian. Starting 
from the energy surface 2?rmf(/3, 7) calculated with the 
DD-PC1 plus separable-pairing functional, each point on 
the (/3, 7) plane is mapped onto the corresponding point 
on the energy surface calculated in the IBM, referred to 
hereafter as Eyqm{P b , 1 b) , using the method proposed 
in Ref. [27} . Here /3b and 75 denote the boson images of 
the quadrupole deformation parameters (3 and 7, respec- 
tively, that are used as constraints in the self-consistent 
RHB calculation and appear as variables in the collective 
Hamiltonian. The boson images /3 b and 75 are related 
to (3 and 7 through the proportionality /3b oc /3, and the 
equality 75 = 7, respectively [2(||27|. This mapping pro- 
cedure is used to determine the strength parameters of 
the IBM Hamiltonian. 

We consider the IBM-2 model [23(: the number of pro- 



ton (neutron) bosons, denoted by {n v ), are assumed 
to equal half the number of valence protons (neutrons). 
In the consistent-Q formalism [22[ the IBM-2 Hamilto- 
nian reads: 

H\bm = e(«d7r + n dv ) + kQk ■ Q v , (4) 

where fidp — dl ■ d p (p — ir or v) and Q p = 
[spdp + d p Sp]^ + Xp[dldp] denote the (i-boson number- 
operator and the quadrupole operator, respectively, e 
and k are coupling constants. The parameters Xn,u in- 
side the quadrupole operators are quite relevant to deter- 
mining whether a nucleus is prolate or oblate deformed. 

The bosonic energy surface £ibm(/3,7) corre- 
sponds to the classical limit of the Hamiltonian 

^ibm: EibmWb,1b) = (*(/?b, 7b)I#ibm|* 7b)>- 
|*(/?b,7b)) denotes the boson coherent state 

0: |*(/3b,7b)) c< U p =,,M + Pp cos ^p4o + 
-^f3pSm r y p (d P+2 + dj 3 _ 2 )] rlf ' |0) , up to a normalization 
constant. Here |0) is the boson vacuum, and the vari- 
ables /3 p and 7 P are the corresponding polar deformation 
parameters. As in our previous studies [HI, it 
is assumed that /3„ = /3„ = /3b and 7^ = 7„ = 75. 
The analytical form of -Eibm(/3b>7b) can be found in 
Refs. |2g,|27]. Hereafter we denote the bosonic energy 
surface as _Eibm(/3,7), omitting the indices of /3b and 
7b- 

The boson Hamiltonian -ffiBM, parametrized by the 
microscopically calculated coupling constants, is diago- 
nalized in the M = boson space. Here M denotes the 
^-component of the total boson angular momentum L. 
Reduced quadrupole transition probabilities _B(E2) are 
calculated for transitions between the eigenstates of the 
IBM Hamiltonian. 

Here we point out again that the total boson energy 
E\bm{P,i) has been related to the microscopic EDF en- 
ergy surface (total energy). However, for the IBM Hamil- 
tonian -ffiBM one cannot make a distinction between the 
kinetic and potential terms, as in the corresponding col- 
lective Hamiltonian H co \\. Nevertheless, the effects rele- 
vant to both vibrational and rotational kinetic energies 
are assumed to be incorporated into the IBM approach by 
adjusting Ejbm(0, 7) to be as close as possible to the mi- 
croscopic surface -Erhb(/3, 7)- This prescription turned 
out to be valid for vibrational and 7-soft nuclei at mod- 
erate quadrupole deformation [2H |27|. similarly to the 
conventional mapping method of Ref. [23j . For rotational 
nuclei with large quadrupole deformation, however, the 
overall scale of the IBM rotational spectra differs from 
the experimental one [26l . |27[ ■ The discrepancy partially 
arises because nuclear rotational properties, character- 
ized by the overlap of the intrinsic state and the ro- 
tated one, differs from the rotational characteristics of 
the corresponding boson system [47l |. This problem may 
be cured by the recently proposed prescription [47j . in 
which the rotational response (i.e., cranking) of boson 
system is related to the rotational response of nucleon 
system. This procedure goes beyond the simple analysis 
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of the zero-frequency energy surface. In order that the 
boson rotational response becomes equal to the fermion 
(nucleon) response, an additional kinetic term L ■ L has 
to be included in the boson Hamiltonian, with a coupling 
constant determined microscopically [47| . The term L ■ L 
directly influences the moment of inertia of rotational 
band with the eigenvalue L(L + 1). However, the above- 
mentioned problem, concerning the IBM rotational spec- 
tra, does not occur in the considered Pt nuclei, and thus 
one does not need to include the L ■ L term in the present 
case. 

Similar problems with the overall scale of the rotational 
spectra are also encountered in the collective Hamilto- 
nian model, when the IB formula is used to calculate the 
moments of inertia [1, 0, Hl| . The inclusion of an addi- 
tional scale parameter is often necessary because of the 
well known fact that the IB formula predicts effective 
moments of inertia that are considerably smaller than 
empirical values. More realistic values are only obtained 
if one uses the Thouless-Valatin (TV) formula, but this 
procedure is computationally much more demanding. In 
the present case we have used the IB moments of inertia 
in the calculation of excitation spectra of Pt nuclei, and 
the agreement with experiment is such that no renormal- 
ization of the effective moments of inertia is required. 
This result allows for a direct comparison of the IBM 
spectra to the solutions of the collective Hamiltonian. 



III. BINDING ENERGY SURFACES IN THE 
f3 - 7 PLANE 

Most deformed nuclei display axially-symmetric pro- 
late ground-state shapes, but few areas of the nuclide 
chart are characterized by the occurrence of non-axial 
shapes. One example is the A « 190 mass region, where 
both prolate to oblate shape transitions, and even triaxial 
ground-state shapes have been predicted. 

The left-hand side of Fig. Q] shows the self-consistent 
RHB quadrupole binding energy maps of the 192 . 194 > 19 6p^ 
isotopes in the (3 — 7 plane, calculated with the DD-PC1 
energy density functional. The energy surfaces are 7-soft, 
with shallow minima at 7 w 30°. In general the equilib- 
rium deformation decreases with mass number and, pro- 
ceeding to even heavier isotopes, one finds that the en- 
ergy map of 198 Pt has also a non-axial minimum, whereas 
200 Pt displays a slightly oblate minimum jl5|, signaling 
the shell-closure at the neutron number TV = 126. On 
the right-hand side of Fig. [TJ we plot the corresponding 
IBM energy surfaces -Eibm^t), obtained by mapping 
each point of surface -Erhb(/?,7) onto the energy sur- 
face calculated in the IBM, following the procedure of 
Ref. [27[. To be able to compare the low-energy spec- 
tra in the two models, the IBM surfaces are mapped in 
such a way to reproduce the RHB energy surfaces up to 
w 2 MeV above the mean-field minimum. This means 
that the maps shown in Fig. [I] can only be compared for 
values of f3 not very different from the minimum /3 m i n . 



For larger values of /3, that is, for higher excitation en- 
ergies the topology of the RHB surfaces is determined 
by single- nucleon configurations that are not included 
in the model space (valence space) from which the IBM 
bosons are constructed. For large /3-deformations, there- 
fore, one should not try to map the microscopic energy 
surfaces onto the IBM. This is the reason why the IBM 
energy surfaces are by construction always rather flat in 
the region f3 ^> j3 m - m . In the vicinity of the minima the 
curvatures of the IBM energy maps are rather similar to 
those of the original RHB surfaces both in j3 and 7 direc- 
tions. The derived values for the Xtt and Xv parameters 
in Eq.(|4]) satisfy x-n + Xv ~ 0, characteristic for a 7-soft 
energy surface. 

One might notice that the IBM energy maps reproduce 
the value of (3 at the minima predicted by the RHB cal- 
culation, whereas the mapping does not reproduce the 
shallow triaxial minima of the RHB surfaces. The min- 
ima of the IBM maps are either oblate or prolate. This is 
because the IBM Hamiltonian of Eq. ((4]) is too restricted 
to produce a triaxial minimum. In the analytical expres- 
sion for £'ibm(/3, 7) the 7-dependent term is proportional 
t° (Xtt + Xv) cos 37, and this places the minimum either 
on the prolate or oblate side according to the sign of 
(x-k +Xu)- The Pt nuclei considered here do not display 
any rapid structural change but remain 7-soft. This fea- 
ture appears to be independent of the choice of the EDF. 
A recent microscopic calculation using the Gogny-DlS 
EDF [l3| also yielded shallow triaxial shapes, rather flat 
in the oblate region [37|, but quantitatively consistent 
with the present analysis. A similar trend was reported 
in other EDF-based studies of ground-state shapes of Pt 
isotopes |39|, |48, 49] . In the present calculation the RHB 
surfaces become softer in 7 with increasing neutron num- 
ber, and the softest nucleus is 196 Pt. The corresponding 
IBM energy surfaces follow this evolution, but do not 
reproduce the triaxial minima because of the reasons ex- 
plained above. The recent Gogny-EDF calculation [37j 
predicts 192 Pt to be the softest Pt isotope in this mass 
region. 



IV. SPECTROSCOPIC PROPERTIES 

In Fig. [2] we display the corresponding low-energy col- 
lective spectra of 192 > 194 > 196 pt obtained from the collec- 
tive Hamiltonian (middle panels), and the IBM Hamilto- 
nian (panels on the right). The calculated ground-state 
and (quasi) 7- vibration bands are compared to the cor- 
responding sequences of experimental states [501 ] . The 
eigenstates of the collective Hamiltonian in Eq. (fTJ) are 
completely determined by the DD-PC1 energy density 
functional plus a separable pairing interaction, and the 
transition probabilities are calculated in the full config- 
uration space using the bare value of the proton charge. 
Since Hibm in Eq. (j4]) acts only in the boson valence 
space, to calculate the B(E2) values one needs two addi- 
tional parameters: the proton-boson and neutron-boson 
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FIG. 1. (Color online) Self-consistent binding-energy maps of 192 . 194 > 196 pt in the f3 — 7 plane (0° < 7 < 60°), calculated 
with the RHB model using the DD-PC1 functional (left panels), and the corresponding mapped energy surface of the IBM, 
Eibm(/3b, 7b). The IBM total energies are depicted in terms of ft and 7, where /3 oc /3s and 7 = jb (see text for definition). 



effective charges. For simplicity, here we take these ef- 
fective charges to be equal, and in each nucleus nor- 
malize the B(E2) values obtained in the IBM to repro- 
duce the transition probability £>(E2; 2^ — > 0^~) calcu- 
lated with the collective Hamiltonian. Thus we can only 
compare the ratios of the IBM B(E2) values, divided by 
B(E2;2± —> Of), to those predicted by the collective 
Hamiltonian based on DD-PC1, and to available data. 

For the ground-state band, both the collective Hamil- 
tonian and the IBM predict excitation spectra in close 
agreement with experiment. For 192 Pt, in particular, the 
calculated ground-state bands seem to indicate a some- 



what larger deformation than observed experimentally. 
In fact, the theoretical energy ratio R4/2 — E(A~l) / E{2~1) 
is 2.59 with collective Hamiltonian, and is 2.69 with the 
IBM Hamiltonian, compared to the experimental value 
R4/2 = 2.48. A similar trend is also found for the other 
two nuclei. A more pronounced difference between the 
predictions of the two models is found in the E2 decay 
pattern of the ground-state band, particularly in 194 Pt 
nucleus for which data are available up to angular mo- 
mentum 10 + . For the spectrum calculated with the col- 
lective model, the E2 transition rates from the state with 
angular momentum L (L 2) to the one with L — 2 keep 
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FIG. 2. (Color online) Low-lying collective spectra of 
192, 194, 196 p^. nuc i e j ; calculated with the collective Hamiltonian 
based on the DD-PC1 functional and the corresponding IBM 
Hamiltonian, in comparison with available data. For each nu- 
cleus, the _B(E2) values (in Weisskopf units) obtained in the 
IBM are normalized to the -B(E2; 2f — > Of) predicted by the 
collective Hamiltonian. The experimental excitation spectra 
and B(E2) values are from Refs. [Hoi ] and [5l|. respectively. 



increasing as function of L, even though the correspond- 
ing experimental B(E2) values in 192 > 194 pt decrease start- 
ing from L = 6. The trend of the B(E2) values calculated 
with the IBM, on the other hand, is much closer to ex- 
periment. The B(E2)'s decrease in the IBM because the 
model space is built from valence nucleons only, and the 
wave functions of higher angular-momentum states cor- 
respond to simple configurations of fully aligned e?-bosons 
[22J , whereas there is no limit on the angular momentum 
of eigenstates of the collective Hamiltonian. 



A more significant difference between the spectroscopic 
properties predicted by the collective Hamiltonian and 
the IBM is found in the sequence of levels built on the 
state 2f - the (quasi) 7-band. The IBM spectra display 
a staggering of excitation energies above 2+, with the 
formation of doublets (3+ 4+), (5+ 6+), ... etc, whereas 
the collective Hamiltonian yields a regular excitation pat- 
tern consistent with the experimental band. To be more 
precise, the IBM spectra correspond to 7-unstable nuclei, 
and are close to the limit of 0(6) dynamical symmetry in 
which eigenstates of a boson Hamiltonian with the same 
r quantum number are degenerate (40j | . On the other 
hand, the 7-bands predicted by the collective model, as 
well as the experimental sequence, seem to be closer to 
rigid triaxiality J52|. The difference between the collec- 
tive Hamiltonian and the IBM arises probably because 
the shallow triaxial minima of the RHB energy surfaces 
are not reproduced by the mapping onto the IBM total 
energy (cf. Fig. [TJ. The agreement of the IBM (quasi) 
7-band with experiment could be improved by introduc- 
ing additional interaction terms in the IBM Hamiltonian, 
i.e., three-body terms (the so-called cubic terms) [53ll54|. 
Terms of this type will have to be included for a more 
precise analysis and comparison of states above the yrast 
with experimental results. 

A nice feature of the present calculation, particularly 
the one with the IBM Hamiltonian, is that the predicted 
B(E2) values for the transition 2^ — > 2^ are comparable 
to or even larger than those corresponding to 4^" — > 
This result is consistent with the experimental trend, 
whereas in the recent Gogny-based EDF calculation of 
Rcf. i3, the 2\ -» 2\ transitions were much weaker 
than 4^ — > 2\ . The corresponding Gogny energy sur- 
faces displayed pronounced oblate minima in Ref. [37| . 
unlike the present energy maps shown in Fig. [T] 

Finally, in Figs. [3] and 0] we compare the absolute 
squares of the collective wave functions for the yrast 
states Of, 2 



1 ! 



4f, 



and the band-head of the 7-band 
of 192 Pt, calculated in the two models. These quanti- 
ties are proportional to the probability density distribu- 
tions in the (3 — 7 plane. Figure [3] shows the distribu- 
tions Em=-l l($il*(/3,7))| 2 , where |$&) denotes the 
IBM eigenstate for the state with angular momentum 
L and projection M. The wave functions of the yrast 
states are concentrated along the oblate axis, only for 
the state 4^ the maximum of the absolute square is lo- 
cated at 7 ~ 55°, and somewhat larger deviations from 
pure oblate configurations are found for higher angular 
momenta. For the state 2^, on the other hand, the peak 
appears in the triaxial region (7 ~ 35°), and the distri- 
bution is extended more toward oblate quadrupole defor- 
mations. The rather large overlap of the collective wave 
functions for the states 2f and 2\ explains the partic- 
ularly strong 2f — > 2f transitions in this nucleus, and 
similarly in the other two Pt isotopes considered here. 
The corresponding absolute squares of the eigenstates of 
the collective Hamiltonian are shown in Fig. 01 In this 
case already the wave functions of the yrast states re- 
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fleet the 7-softness of the RHB energy surface, and the 
maxima of the absolute squares are found in the triaxial 
region of the j3 — 7 plane. 



V. CONCLUSIONS AND OUTLOOK 

Structure phenomena related to shape evolution cur- 
rently present a very active research field in low-energy 
nuclear physics. Radioactive-beam facilities continue to 
provide interesting new data on shapes in regions of ex- 
otic nuclei far from stability. The variation of ground- 
state shapes is, of course, governed by the evolution of the 
underlying shell structure of single- nucleon orbitals. It is, 
therefore, important to develop microscopic approaches 
that can be employed in quantitative analyses of shape 
phenomena and the resulting complex excitation spectra 
and decay patterns across the entire chart of nuclides. 
Such a framework is provided by nuclear energy density 
functionals (EDFs). 

The advantages of EDFs are evident already at the 
basic level of implementation: an intuitive interpreta- 
tion of self-consistent mean-field results in terms of in- 
trinsic shapes and single-particle states, and the univer- 
sality of EDFs that enables their applications to all nuclei 
throughout the periodic chart. The latter is particularly 



important for extrapolations to regions of short-lived nu- 
clei far from stability. When considering spectroscopic 
applications, the framework of EDF must be extended 
beyond the mean-field level to allow for a systematic 
treatment of dynamical effects related to restoration of 
broken symmetries and fluctuations in collective coor- 
dinates. To calculate excitation spectra and transition 
rates, it is necessary to project states with good quantum 
numbers from the mean-field solution, and also take into 
account fluctuations around the mean-field minimum. 

In this work we have compared two well known mod- 
els that explicitly consider quadrupole collective correla- 
tions, both starting from maps of binding energy cal- 
culated with the same microscopic EDF. The first is 
the generalized collective Hamiltonian for quadrupole 
vibrations and rotations. The dynamics of the five- 
dimensional Hamiltonian is governed by the collective 
potential, the three vibrational mass parameters, and 
three moments of inertia for rotations around the princi- 
pal axes. These functions of the quadrupole deformation 
parameters are determined by constrained mean-field cal- 
culations using a given nuclear EDF. The diagonalization 
of the resulting Hamiltonian yields excitation energies 
and collective wave functions that can be used to calcu- 
late various observables. Calculations are performed in 
the full model space of occupied states (no distinction 
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FIG. 4. (Color online) Same as described in the caption to Fig. [3] but for the eigenstates of the collective Hamiltonian. 



between core and valence nucleons, no need for effective 
charges). The second model considered in this work is the 
well-known and very successful IBM-2. In this approach 
the configuration space is first restricted to valence nu- 
cleons only, and further mapped to the space of s and d 
bosons. To determine the parameters of the IBM Hamil- 
tonian, the energy surface calculated using a microscopic 
EDF, is mapped onto the corresponding boson energy 
surface under certain approximations. One then proceeds 
to calculate the excitation spectra and wave functions in 
the laboratory frame. To calculate transition probabil- 
ities, however, one needs to adjust the effective boson 
charges. Here this has been done so that in each nucleus 
the calculated B(E2; 2^ — > 0^~) coincides with the value 
obtained using the collective Hamiltonian. 

The two models have been compared here in a study of 
the evolution of non-axial shapes in Pt isotopes. Start- 
ing from the binding energy surfaces of 192 . 194 ^ 196 pt j 
calculated with the DD-PC1 energy density functional 
plus a separable pairing interaction, we have analyzed 
the resulting low-energy collective spectra obtained from 
the collective Hamiltonian, and the corresponding IBM- 
2 Hamiltonian. The calculated ground-state and 7- 
vibration bands have been also compared to the corre- 
sponding sequences of experimental states. Both models 
predict that excitation energies and B(E2) values are in 
agreement with data. In particular, we notice the ex- 



cellent result for the predicted excitation energy of the 
band-head of the 7-band, as well as the good agreement 
with the experimental B(E2) values for transitions be- 
tween the 7-band and the yrast band. 

There are also significant differences in the predictions 
of the two models. With the present form of the IBM 
Hamiltonian, restricted to two-body boson interactions, 
its expectation value in the boson coherent state does 
not reproduce the shallow triaxial minima of the binding 
energy maps predicted by the constrained self-consistent 
mean-field calculation using DD-PC1. Since the mapped 
IBM energy surface is 7-soft rather than triaxial, the re- 
sulting spectra display a staggering of excitation energies 
above 2+, with the formation of doublets (3+ 4+), (5+ 
6+), ... etc, in contrast to the regular excitation pattern 
observed in experiment and reproduced by the collective 
Hamiltonian. This problem could be solved by including 
three-body boson terms in the IBM Hamiltonian, and 
work along this line is already in progress. When con- 
sidering the calculated B(E2) values for transitions in 
the ground-state band, the IBM reproduces the grad- 
ual decrease of transition rates with angular momentum 
for L > 6, reflecting the finiteness of the valence space. 
On the other hand, even though the collective Hamil- 
tonian predicts parameter-free _B(E2) values in excellent 
agreement with experiment for transitions between low- 
spin states, the calculated transition probabilities keep 
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increasing with angular momentum, in contrast to data. 

Both models are based on binding energy surfaces cal- 
culated at zero rotational frequency. In general this leads 
to effective rotational moments of inertia that are lower 
than empirical values, that is, the calculated rotational 
bands are stretched in energy compared to experimental 
bands. In the collective Hamiltonian the moments of in- 
ertia can be improved by including the Thouless-Valatin 
dynamical rearrangement contributions. For the IBM 
Hamiltonian one needs to include the kinetic rotational 
term [471 ] . and perform the mapping of microscopic en- 
ergy surfaces calculated at finite values of the rotational 
frequency. We have already started with the implemen- 



tation of these modifications in our current version of the 
collective Hamiltonian based on relativistic EDF, and in 
the IBM Hamiltonian. The comparison of the improved 
models will be the subject of a future study. 
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